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Abstract 

One-dimensional quasilattices, namely, the geometrical objects to represent 
quasicrystals, are classified into mutual local-derivability (MLD) classes. Be- 
sides the familiar class, there exist an infinite number of new MLD classes, 
and different MLD classes are distinguished by the inflation rules of their 
representatives. It has been found that electronic properties of a new MLD 
class are characterized by the presence of marginal critical states, which are 
considered to be nearly localized states. 
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Quasicrystals (QCs) have aperiodic ordered structures, which are different from either 
of crystalhne and disordered materials. It has been reported that QCs exhibit quite unique 
physical properties which should be closely connected with the structure of QCs. A 
sufficient understanding in this respect is, however, yet to be attained. For instance, in 
order to understand the transport properties of QCs, further systematic studies on the 
electronic properties of various types of QCs are necessary. 

The atomic positions of an idealized QC form a quasiperiodic object called quasilattice 
(QL), the classification of which is the principal subject of the crystallography of QCs. One 
of the important features of QLs is the self-similarity, which is closely connected with the 
physical properties. Thus, in the first part of this letter, we shall develop a new classification 
scheme of QLs based on the self-similarity. One can argue that the different classes of QLs 
correspond to different universality classes with respect to the electronic properties. We next 
show that there exists a new class of QLs whose electronic properties are characterized by 
the presence of marginal critical states, which can be thought of as being virtually localized. 
This feature is essentially different from the case of the standard class of QLs on which most 
of the previous studies have been done. The transport properties of the new class may be 
markedly different from the conventional one. 

We confine ourselves to one-dimensional (ID) cases. A ID QL is obtained by projecting 
a subset of a 2D lattice A, the mother lattice, onto a ID subspace E\\ [0,0. The line E\\ is 
taken to be parallel to rai + sl2 with ai and sl2 being the primitive lattice vectors of A and 
T a quadratic irrational such as |(1 + V^) (the golden mean), 1 + \/2 (the silver mean), 
etc. The subset is taken as A fl S with E being a parallel strip to E\\. If the scales of £"11 
and E±, the physical space and its orthogonal complement, are chosen appropriately, both 
the projections Ay and A_|_ of A onto E\\ and E^, respectively, are given by the dense set 
Z[r] = {p + qr \ p,q E Z}. Thus, a ID QL is a discrete subset of Ay, where its point density 
is proportional to the window W, i.e., the width of the strip S. According to the relative 
position of S to A, an infinite number of QLs are possible with the same W. But they form a 
single local-isomorphism class, and their differences are irrelevant to the physical properties. 

For a generic W, one obtains a ternary ID QL composed of three types of lattice spacings, 
s, m and 1, referring to short, medium and long spacings, respectively; they satisfy |1| = 
|s| + |m|. If W is reduced, one obtains another QL as a subset of the original QL. The new 
spacings, s', m' and 1', are composites of the older ones. For an appropriate choice of the 
reduction of W, the older QL turns out to be a uniform decoration of the newer, i.e., the 
decoration of 1', for example, is common for all I's. The three decorations of the spacings 
form a substitution rule. The newer QL is then called a subquasilattice (SQL) of the older. 

If W is reduced by the factor r~" with n being any positive integer, the resulting QL 
is locally isomorphic with the scaled version of the original QL with the factor r" [Q. 
The original QL is called self-similar if the new QL is an SQL of the original. Here, the 
substitution rule, combining the three spacings of the original QL with the newer ones, 
represents the self-similarity; it is called the inflation rule. It can be shown that a QL can 
be self-similar if and only if its window W belongs to the quadratic field Q[r] = {r+sr | r, s G 
Q} DH; that is, the window must be rational. The inflation rule as well as the minimum 
power n for the ratio of self-similarity, r", depends on a number theoretical property of W. 
The following argument will be focused on this rational case. 

Several QLs associated with the silver mean are presented in Table |. The first two are 



2 



binary but others ternary. The QLs except for "C" and "E" have self-similarities whose 
ratios are presented in the 6th column. The QL "C" is not self-similar but is changed to the 
once-infiated version of "A" by the substitution rule in the 7th column. A similar relation 
holds between the pair {D, E}. The QL "A" is an SQL of "B" with the substitution rule, 
s' = s, 1' = sm, and the QL "F" is an SQL of "G" with s' = m, m' = 1, 1' = msm. Though 
the two QLs "G" and "H" are similar, there exists no inflation rule which combine them. 

We introduce an important relationship between QLs: two QLs are mutually locally- 
derivable if all the sites of one of them are determined locally from the structure of the other 
and vice versa 0. All the QLs derived from a single mother lattice are classified into an 
infinite number of mutual local- derivability (MLD) classes. The MLD class to which a given 
QL belongs is determined by a number theoretical property ofW G Q[t]. We only present a 
necessary condition: if two QLs belong to a single MLD class, the relevant two windows have 
a common denominator when they are represented as simple fractions in the quadratic field. 
A QL and its any SQL belong to the same MLD class. The eight QLs in Table | are divided 
into three MLD classes: {A, B, C}, {D, E}, {F, G, H}. It can be shown that i) every MLD 
class includes at least one self-similar member, which can be taken as a representative of the 
class, and ii) every QL in the MLD class has an SQL which is similar to the representative; 
that is, any non-self-similar QL is a uniform decoration of a self-similar QL. 

The structure factor of a QL is composed of Bragg peaks, whose intensities are determined 
by the size of the window. For example, the three QLs, "C", "D" and "F" in Table |, 
have windows of similar sizes, and their structure factors are not so much different. They 
nevertheless belong to different MLD classes. 

We shall call a QL to be type I or II according as its window belongs to Z[r] or not, 
respectively. All the type I QLs for a given r form a single MLD class, as seen in Table |. 
Previous investigations on the electronic properties of ID QLs have been almost exclusively 
done on the basis of models on type I QLs It has been proved that a type I QL and its 

any decoration belong to a common universality class of electronic properties |TD|] , which are 



dominated by the structures of the inflation rule. On the other hand, we may expect that 
type II QLs belong to different universality classes, because their inflation rules are different 
from that of the type I class. 

We now present a brief review on the electronic properties of type I QLs. The main 
observations, obtained from the case of the Fibonacci lattice as well as some of its associates, 
are 0: i) the energy spectrum is purely singular- continuous or, equivalently, fractal-like, 
and ii) all the eigenf unctions are critical, i.e., neither extended nor localized in the usual 
meaning. The energy spectrum and the eigenfunctions have also been found to exhibit self- 
similar structures, which can be directly related to the inflation rule of the underlying QL 
by a real-space renormalization-group approach . 

The energy spectrum of any homogeneous ID system obeys locally a scaling law. The 
local scaling at the reference energy Ej. is represented by the scaling exponent a = a{Ej.) 
satisfying < a < 1 0. While a represents the local dimension of the energy spectrum, 
it also characterizes the localization character of the eigenfunction of the energy level E^.. 
In particular, an isolated energy level of a localized state has a vanishing exponent, a = 0, 
while inside an absolutely continuous spectrum, which confirms extended states, a = 1. A 
purely singular-continuous spectrum being characterized by fractional exponents confirms 
critical eigenfunctions. The energy spectrum of a type I QL is, in general, a multifractal 
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TT| , and is characterized by the /(a) spectrum. The support of the f{a) is an interval 
[amin,(ymax] with < amin < Oimax < I] a rauges from amin to Umax, SO that every energy 
level is characterized by fractional power-law scaling |]^. 

Although the above observations are common in the case of type I QLs, electronic proper- 
ties of type II QLs have been scarcely investigated. In the following, we show some numerical 
as well as analytical evidence that the type II QL listed as "D" in Table H and shown in 
Fig. |I| indeed exhibits a new scaling property. Let us take a binary atomic chain, which 
is obtained by decorating this QL as follows: i) atoms of type X are located on all the 
lattice points, ii) one atom of type Y is located on each spacing of type m, and iii) a pair 
of type Y atoms is located on each spacing of type 1. The densities of the two types of 
atoms in the chain are equal. More remarkably, the chain is invariant against the exchange 
of two types of atoms, X and Y. We employ the tight-binding model on this atomic chain, 

j-i + j + j+i = ji where we assume that t = —1 for the transfer integrals and 
Vx = and Vy = \^ for the relevant site energies. 

Due to the aforementioned symmetry of the atomic chain, the energy spectrum becomes 
symmetrical which is confirmed by the numerical result as shown in Fig. ^ The energy 
spectrum is divided into two equivalent clusters, and the left cluster exhibits a hierarchical 
trifurcation, which will reflect the nature of the leftmost level, i.e., the ground state. The 
trifurcation can be understood in terms of a perturbational real-space renormalization-group 
approach similar to that introduced by Niu and Nori for the case of the Fibonacci lattice 
P). The approach is based on the recursive structure of the QL and assumes V ^ 1. In 
the zeroth approximation, the spectrum consists of two equally degenerate energy levels at 

= and V , while the eigenstates are the Wannier states themselves. To see the splitting 
of the level originating from X atoms, all the Y atoms are decimated as shown in Fig. |^. 
This yields three types of effective transfer integrals between X atoms, = —1, tm ~ —^/V, 
and t\ ^ —1/V'^, which are all negative and satisfy the inequalities: 

\ts\ > |tm| > \tl\. (1) 

In the first approximation, we take only the leading effective transfer integral into 
account. Since the type s spacings are isolated, the energy spectrum consists of three sub- 
levels a.t E = and E = ±ts (= =Fl); the central level is derived from isolated atoms, 
while the two satellites from isolated diatomic "molecules". The weights of the three sub- 
levels have the ratios which agree with the numerical results. To see the 
trifurcation of the level derived from the molecular bonding states, we shall consider the 
"molecules" to be the "atoms" of the second generation. Due to the self-similarity of the 
QL, there appear three types of effective transfer integrals between the new "atoms" as 
shown in Fig. ^: 

t>ltm, 4~^ii, andt;^^|, (2) 

where the numerical factor 1/2 appears as the square of the bonding amplitude l/V^. Since 
the new parameters satisfy the same inequalities as Eq. (|l|), we can return now to the 
point after Eq. (|I|) and continue the discussion recursively. Thus, we obtain a hierarchical 
trifurcation spectrum. The ground state is understood to be a hierarchical composite of 
molecular bonding states, as shown in Fig. 11. 
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To make a quantitative argument, we introduce two positive parameters as the ratios 
between the effective transfer integrals of the nth generation: 

t-s tin 

Eq. d^) yields the two-dimensional map, fn+i ~ Qn, dn+i ~ /n5'n, which converges to 
the origin of the /(/-plane. This means that the "atoms" are increasingly isolated as n is 
increased. We can linearize this 2D map with new variables, Infn and In^f^, and obtain the 
asymptotic expression 

fn ~ gn-1 ~ exp [-CT^] = CXp [-Ct'"'], (4) 

with tg = (1 + a/5)/2, i.e., the golden mean, and u = Inrc/ Inr ^ 0.546. Note that tq is the 
leading eigenvalue of the 2x2 matrix associated with the linearized map and c a positive 
constant depending on V. 

The band width Wn of the subcluster of the nth generation is estimated to be w„ ~ 2|)f:^")|. 
The relation ^ ^t^^ as in Eq. (||) can be used to prove Wn+i/Wn ~ i(:^"+^V4"^ ~ |/n ~ 

exp [— cr"'^], which tends to zero as n — oo. This is observed in Fig. ^ Then, it can be 
proved that Wn ~ exp [— c'r"^] with c' being another positive constant. Since the size L of 
the "atom" of the nth generation is proportional to r", we find w„ ~ exp [— (L/^)^], i.e., a 
stretched exponential, with ^ being a characteristic length. This is manifestly different from 
the case of usual critical states exhibiting a power-law, Wn ~ L^^/", and has never been 
observed in energy spectra of type I QLs. Note that a power-law scaling could be observed 
if the ratio Wn+i/wn had tended to a finite value. 

The above argument proves that the exponent a of the ground state vanishes. This does 
not mean, however, that the ground state is localized, because the ground state level is not 
isolated from other levels. The length ^, which decreases as V is increased and diverges 
as y — i> 0, represents the size of almost isolated "atoms". The presence of a characteristic 
length implies that the ground state wave function, shown in Fig. ^, is not self-similar, 
although it still remains critical. We may call the ground state a marginal critical state, 
presence of which is closely connected with the self-similarity of the present QL. 

The exact real-space renormalization-group formalism proposed by Ashraff and Stinch- 
combe can also be adapted to study the ground state properties. There appears a 7D map 
but its asymptotic behavior is dominated by its behavior in a 2D subspace corresponding 
to the /(/-plane. The map being reduced to the subspace is essentially the same as above. 
It has also been numerically confirmed that each energy level at a band edge, where energy 
levels only accumulate from one side, gives the equivalent fixed point of the 7D map to that 
of the the ground state level. 

The presence of marginal critical states causes vanishing of the left end amin of the sup- 
port [amin, Cimax] of the relevant /(a) spectrum. This is an important feature of electronic 
properties of type II QLs, and has also been confirmed for several other type II QLs in- 
cluding the one whose ratio of self-similarity is Tq. The exponent u introduced above is a 
proper number to the relevant MLD class but satisfies < u < 1. It is surprising that QLs 
derived from a single mother lattice can belong to different universality classes of electronic 
properties. 
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The binary sequence of atoms in the atomic chain used to investigate electronic properties 
of type II QLs is found to be one of the circle sequences investigated in detail by Aubry et 
al Some of other type II QLs can also be transformed to circle sequences if they are 



decorated appropriately, but a general relationship between the type II QLs and the circle 
sequences awaits a further investigation. We should remark, however, that an indication of 
an unusual electronic property of a circle sequence was first reported by Luck ||T3[| . 

The present report is summarized as follows: i) ID QLs can be classified into an infinite 
number of MLD classes, which are distinguished from each other by the inflation rules of 
their representatives, and ii) there exists a "type 11" MLD class whose electronic properties 
are characterized by the presence of marginal critical states. 

A full account of the exact real-space renormalization-group treatment of marginal crit- 
ical states will be presented elsewhere. The present results can also be extended to 2D and 
3D QLs because these QLs are related to ID QLs through the Ammann bars or planes 0, 
and the results will be reported elsewhere. 
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FIGURES 



FIG. 1. A ternary ID QL given as "D" in Table I. It is self-similar with the inflation rule in 

the table. The ratio of the self-similarity is the silver mean, t = 1 + The frequencies of the 
three spacings s, m, and 1 are shown to be proportional to : : (2t~^ -|- = 1). 



FIG. 2. A singular-continuous energy spectrum of the type II QL given in Fig. 1. It was 
obtained numerically from a finite approximant composed of 1970 atoms with V = 5.0. The left 
cluster of energy levels exhibits trifurcating behavior when it is expanded successively. 



FIG. 3. A perturbational real-space renormalization-group treatment of the ground state. 
"Molecules" are separated by three types of spacings, which are arranged similarly to the un- 
derlying QL. 



FIG. 4. The ground state wave function of the approximant composed of 1970 atoms with 
V = 1.0. is the probability amplitude on the jth site. This figure exhibits the molecular 
bonding state of the eighth generation, where hierarchical internal structure is observed. 
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TABLES 





W 


|s| 


|m| 


|1| 


ratio 


inflation rule I 


[01 substitution rule) 


A 


V2 


1 


— 


r 


r 


s' = 1, 1' = sll 




B 


I + V2 


1 






r 


s' = sm, m' = 


sms 


C 


3- ^2 


1 


V2 


r 




( s' = sm, m' 


= sll) 


D 


l + ^\/2 


1 


^/2 


r 


T 


s' = sm, m' = 


si, 1' = sll 


E 




1 




r 




( s' = sm, m' 


= ssm, 1' = slsm ) 


F 


\ + V2 


1 


V2 


T 




s' = smsl, m' 


= smsmsl, 1' = smsUsmsl 


G 


l + ^/2 


r-i 


1 






s' = msm, m' 


= mlmmsm, I' = mimlmmsm 


H 




1 


r 


V2r 


r2 


s' = msm, m' 


= mlmmsm, 1' = mimlmmsm 



TABLE I. Several ID QLs associated with the silver mean, r = 1 + \/2- The first three belong 
to the type I MLD class, while others to type II MLD classes. 
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